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Appendix A 

Derivation of the Frequentist and Bayes-modified Likelihood Functions 

The notation used in this Appendix is the same as what has been introduced in the 
main portion of this application. In what follows, we make the following assumptions: 

(i) The visibility of a particular pilot at (x,y) g ^4 at a given time i is independent of 
whether the pilot was visible at any time 5</, and is independent of / as well. 

(ii) The visibility of pilot / at time i is independent of the visibility of pilot jQ^i) at 

time /. 

(iii) Throughout the duration of the prediction process, the mobile is relatively 
stationary, i.e., the mobile unit remains in its initial sub-cell. 

Defining L^(x,y)-l ,V (x,y) gA, a recursive form for the exact likelihood 
function through the first s measurement epochs for the unknovm location of the mobile 
based on (i)-("0 is 

U(x,y) oc L^''ix^y)Yl[0,(x,y)rHl ' 0,(x.y)t'^' , (x,y) e A. (Al) 

The functions Oy(x,y) in (Al) are unknown, so the exact likelihood function is not 

calculable. Replacing the O^jix.y) with their approximations ^^(x,^), defined by equation 

(B9) in Appendix B, gives the definition of the frequentist likelihood function L\^(x,y) 
shown in expression (1). 

The frequentist likelihood fiinction (1) was derived under the assumption that a 

family of functions Oij {x,y) accurately approximate the unknown functions 0ij {x,y) . The 

approximations ff^j^x^y) are based on the RF model described in Appendix B. The RF 
model is an approximation to a very complex stochastic process. Bayesian methods allow 
the uncertainty associated with 0^j{x,y) to be incorporated into the model framework. 

Rather than assuming 5y (x,>^) = (x,^) , we can model the parameters 6^j(x,y) with 
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independent prior distributions that have mean values equal to 0ijix,y). For each 
(x,y)GA, we use a prior distribution for 0yix,y) that is a beta distribution with 
probability density function given by 

5 where ^v] is the complete beta function. Note that the two parameters, a,jix,y) and 
fiy{x,y) , to be determined in what follows, depend on both the possible location {x,y) e A 
and the pilot ij ^K. Whereas the frequentist likelihood is defined by implicitly assuming 
Q ^AGfj{x,y) = e^ix,y)'\ = \, the Bayes-modified likelihood is defined by taking the 

y expected value of the exact likelihood (Al) with respect to the beta distributions defined by 

10 (A2). The recursive form of the exact likelihood (Al) does not amend itself to evaluating 
the expected value as easily as the non-recursive form of the exact likelihood which is 
m L\x,y)o,Yl{e,j(x,y)-f'\\-e,ix,y)r'^ . {x,y-)^A (A3) 

where is the number of times pilot ij was visible through the first s measurement epochs. 
Evaluating the expected value of (A3), with respect to the beta distributions defined by (A2) 
15 gives the Bayes-modified likelihood fixnction 

L's^ix^y) E\Y[[e,j{x.y)f\.i-e,jix,y)rA 

\lJ^K J 

20 



1=5 



It is easily verified that equation (A4) has the recursive form shown in the functional 
expression (4). 

25 The definition of the Bayes-modified likelihood function will be complete once 

values for a^jix.y) and Pij{x,y) are specified. Determining what values to use for 
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<^ij i^yy) fi^ix^y) is an example of the classic Bayesian dilemma - how to specify a 

prior distribution? An underiying principle in specifying a prior distribution is that any 
available information about the parameter of interest, in this case 6ij{x,y), should be 

reflected in the prior distribution that is used. In this case, we have an approximation of 

O^jix^y) from the RF model, namely 6yix,y), which we will use as the mean of the prior 

distribution for 0.j(x,y) , Accordingly, we have the following constraint: 

— - — ' ^ . T = ^ijix^y). (A5) 

cc,j(x,y)-^fi,j(x,y) 

Since we have two unknovms, (x, j^) and Pij{x,y), we need one additional constraint to 
uniquely determine both a^jix.y) and fiij{x,y). To obtain the second constraint, we 

impose that the prior be as vague as possible, subject to the first constraint. To implement 
this second constraint, we will maximize the variance of the prior, subject to the fixed mean 
value. The second constraint is consistent with the quantity of the prior knowledge we have 
about 6ij{x,y) and except for anticipating a value near 9ij{x,y), we know nothing else 

about its potential value. We also restrict attention to beta distributions which have 
bounded density functions, implying ay{x,y)>\ and p^j{x,y)>\. The restriction to 

bounded density functions only rules out "U-shaped and "Slide-shaped" beta distributions 
which would not typically be intuitive representations of the prior information on 0fj{x,y) . 

The problem of selecting a^jix.y) and p^jix.y) has now been sufficiently constrained so 

that it can be solved by using the following Lemma. 
Lemma 

Consider the family of beta distribution with parameters a ^ 1 and b>\. The values of a 
and b which maximize the variance, subject to the distribution having a given mean value 
//,are: a = l and A = , if /i<l/2; a = /i/(l-//) and 6 = 1, if /i>l/2. 

Proof of Lemma 

The variance of the beta distribution is given by = ah l[{a + b + V){a + by\ Since we 
have // = a/(a-i-6), we can write b = a{\- jj)/ fu, and hence = ^^(l- ^)/(a-\-/i). To 
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maximize the variance, we need to minimize a, subject to the constraint that both a > 1 and 
b>\. Clearly, if //<l/2 then we can choose a = l and 6 = (1 -//)///> 1 will be 
ensured. Likewise, if //> 1/2 then the smallest we can choose a and still have 6> 1 is 
a = /i/(l-//)^l. 

Applying the Lemma v^th a = ay(x,y), b = /}.j{x,y) and ju = 0ij{x,y) gives the 
formulas in equations (2) and (3) for a^jix.y) and Pijix.y) . 



